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Abstract. A linear polyomial non-negative on the non-negativity domain of 
finitely many linear polynomials can be expressed as their non-negative linear 
combination. Recently, under several additional assumptions, Helton, Klep, 
and McCuUough extended this result to matrix polynomials. The aim of this 
paper is to study which of these additional assumptions are really necessary. 



1. Introduction 

We are interested in matrix generalizations of the following variant of the Farkas 
lemma. 

Theorem 1. Let /i, /2, ■ • • , A- be linear polynomials in n variables, i.e., fi{xi, X2, 
. . . , Xn) = + a'^i^xi + . . . + an^Xn, where a^*'' e R for i = 1,2, ... ,k, j = 
0,1,..., n. Let K = {x eK" \ f,{x) >0 ior alii = 1,2,. ..,k}. If f is another 
linear polynomial in n variables, for which /|_fs: > holds, then there exist non- 
negative constants Ci, .such that 

f = Cq+ Cl/l + C2/2 + . . . + Cfc/fc. 

We write M''^'*[a;] (resp. §R'^^''[x]) for the set of all polynomials whose coeffi- 
cients are d x d (resp. symmetric d x d) matrices. The evaluation of a linear poly- 
nomial L{x) = Po + T,i=i P^^^ G SR''^''[a;] at X = {Xi,X2, . . . ,X„) e (SM"^™)" 
is defined as 

n 
1=1 

with the usual tensor product of matrices. Let 

DLim) = {Xe (§M"^")" I L{X) h 0}, 

00 

Dl ^ \J DL{m). 

■m—l 

The following generalization of Theorem [T] was obtained by Helton, Klep, and 
McCuUough in [2 . It is a special case of their Theorem 6.1. 

Theorem 2. Suppose Li = ^ + Pi^i ^ §R'^^''[x] is a manic linear polynomial 

and the set Dl-^{1) = {x G R" | Li{x) >^ 0} is bounded. Then for every linear 
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polynomial L2 = Bo + X^ILi Ri^i £ SK^^''[a;] with ^211)^^(1) >- 0, there are matrix 
polynomials Aj e M^^^ [a;] and Bk £ W^^^ [x] satisfying 

L2 = J2A*A,+J2BlL^Bk. 

3 k 

Note that this result also covers the case of several constraints; simply take Li 
to be their direct sum. 

It was already shown in [2^ that the matrix polynomials Aj , Bk need not be 
constant unless the condition L2\dl-^{i) ^ is replaced with L2\dl-^ ^ 0; cf. [2, 
LP-satz]. 

The aim of this paper is to study the necessity of the following assumptions in 
Theorem [2| 

(1) Boundedness of Dl-^{1): By Example 1 this assumption cannot be re- 
moved. 

(2) Monicity of Li: In Theorem [3] we prove that for diagonal Li monicity 
can be removed from Theorem [2l The general case remains open. 

(3) Strict positivity of L2\dl-^{i)'- In Section 4 we show that this assumption 
can be replaced with L2\dl-^{i) h in the following special cases: 

• In the one- variable case (even if £)ij(l) is unbounded). 

• When the span of the coefficients of Li is closed under multiplication. 
The general case remains open. However, by Example 2 simultaneous gen- 
eralization to non-monic Li and non-strict L2\di^^{i) is not possible. 



2. Boundedness of -Dli(I) 

The assumption of boundedness in Theorem [2] cannot be removed, because of 
the following example: 

Example 1. For linear polynomials 



Li 



l + xi 





1 





Xl 





X2 






1 + 2:2 



L2 



4 



3 

- 5x2 



we have that the set Z)li(1) is contained in the set D^^i^) '■— {(2:1, 2:2) £ | 
^2(2:1, 2:2) >- 0}, but L2 cannot be expressed as ^'j^j + J2k ^k^'^^Bk- This im- 
plies that the boundedness of Dl-^{1) is needed in Theorem 2. 

Proof We have DL^il) = {(2;i,X2) G | 2;i > -l,xi-|-2;2 > -1,2:2 > -l}- On 
the other hand, the fact 1 + ^xi > on Dl-^{1) together with 



det(L2) 



1+3.1 



1 + 772:2 - - = 



2 1, 

3 + 3(1 + ^1) 



2 1, 

3 + 3(1 + ^2) 



2 1 
+ - (1 + 2:1 + 1 + 2:2) + - (1 + 2:1) (1 + 2:2) 



^ + ^ (1 + 2^1 + .2) + ^ (1 + 2:1) (1 + 2:2) , 



A NOTE ON A MATRIX VERSION OF THE FARKAS LEMMA 



3 



where the second and the third summand in the last hnc arc non-negative on 
-DliCI), gives L2\dl^{i) >- 0. 

Now we are going to show, that if L2 can be expressed as L2 = + 
B'^LiBk, then Aj and can be assumed to be constant matrices. Let us 
denote Aj , B^ as 



P^'\xuX2) R['\xi,X2) 



P^'\xi,X2) lii>{xi,X2) 



U), 



Pi\xi,X2) ri"'^(xi,X2) 

P2\xi,X2) rf'{xi,X2) 



pi^\xi,X2) rf'{xi,X2) 



Comparing the entry (11) in A*jAj + B^LiB). and L2 gives 

^ [{Pi'\x^,X2)r + {PPix^,X2)r) + J2 {ip[''\x^,X2)r{l + X,) 



+ 



+ ip^^\xi,X2)fil + XI+ X2) + ipi^\xuX2))\l + X2)) ■ 



•1 



1 



-Xi. 



By observing the monomial of the form Kx'!^, n€N,K^Q,i = l,2, oi the highest 
degree on the left side, it can be seen, that monomials Ax^,n € N, A 7^ 0, i = 1,2, 



do not appear in any P^^^ nor p'f . With the same reasoning applied to the entry 



(22) not even in and rf . Further on, since the monomial Kx^x^, K ^ 0, 

m, n G N, from P^ ' ov pl , does not multiply into monomial K or Kxi, if ^ 0, it 
is not needed to satisfy the upper equality. Similar reasoning can be applied to the 
other entries, hence WLOG Aj, Bk arc constant matrices. 

The comparison of coefficients in ^^^j + X^fe Bk^iBk and L2 in a constant- 
matrix case gives the following equalities: 

Entry (11): 



(1) 



X2 : 



E 

k 



{{P: 



i'^r + ipi'^r 







(fe) 
P2 



(fe) 

■Pa 



0, Vfc 



(2) 



Xi : 



k 



j \i=l J k 

Entry (12)(=entry (21)): 
(4) X2: E(P2'=Vf^+^'i'=^r-f)=0 



(5) 



X,: Y,(pfrf^pfrf)=^ 



(6) 



1 : 



E E^. 
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Entry (22): 



(7) 
(8) 



k 



0, V/ 



(9) 



1=1 



We will see, that the upper equalities cannot be simultaneously satisfied. 

= 0. We use this in [B] and get 



(fc) (fc) 



From [T] and [5] we conclude J2k 

(E -=1 ^i'^^F^) - !■ Using [2 and [D in [3] gives 

Similarly using [7] and [S] in El gives J2j (ELi(^P^)^ 
The following chain of (in)equalities should hold: 



2,2 
3 3 



EE- 

j 1=1 



>EEN'^^Pl^EE^^^^i^^ = i' 



j 1=1 



i '=1 



where the first inequality follows from AG-inequalites applied to pairs |(-P/"''')^, 



2 

We conclude § > |, which is obviously a contradiction. □ 



3. MONICITY OF Li 

In this section we show, that for diagonal Li monicity in Theorem [5] can be 
removed. In Proposition[T]we first prove the case of the set Z?Lj (1) being a singleton 
and then also the other cases of bounded 13^^(1) in Theorem[31 

Proposition 1. Suppose Li G §M'^^''[a::] is a diagonal linear polynomial and the set 
= {a}. Then for every linear polynomial L2 G §R^^^[a;] with L2\dl-^{i) ^ O7 
there are matrix polynomials Aj G R^^^ [a;] and Bk G R'*^^ [a;] satisfying 

J k 

In the proof we will use the following proposition: 
Proposition 2. For every A G SE^^^ there exist Bk G §R^^^, such that 



X 

-X 



Bk = Ax. 



Proof. According to the well-known fact every real symmetric matrix is real congru- 
ent to a diagonal D with elements 1,-1 and on the diagonal, i.e., A = J^k B^DBk, 



where £>, Bk G SI 



Dx can be constructed from 



X 

-X 



with the aim of 
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equalities: 



Ei', 



-Ei, 



[ e. 
[ e, 



" 1 







e* 





-1 







" 1 













-1 




e* 



where e, denotes the standard . 



vector. 



□ 



Proof of Proposition]^ Up to translation we may assume 13^^(1) — {0}. For a 
polynomial 



Xi 







-Xl 









-Xr. 



we have Z?^^ (1) = {0}. After applying Theorem [T] to diagonal entries of Li and each 
diagonal entry of Li, it follows Li = A*Aj B^LiBk for some constant Aj, 
Bk- Therefore it suffices to find Aj,Bk, such that L2 = + Efc BlLiBk. 

If we write L2{x) — Rq + J27=i ^i^i^ then L2(0) ~ Ro h 0- So there exists 
A, such that Rq — A* A. According to Proposition [5J RiXi can be expressed in a 
Xi 



desired way with 



—X, 

Since i was arbitrary, we are done 



hence also with Li. 



□ 



Now we will extend Proposition [T] to the general case. One additional lemma 
will be needed for that. 

Lemma 1. Suppose L = Pq + X^iLi Pi^i £ §K''^''[a;] is a linear polynomial and 
let be an interior point of the set Dl-^{1). Then there exists a monic linear 
polynomial L = I + X]"=i Pi^i ^ SM''^'^[x], where d < d and -Dl(I) = £'^(1), such 
that L = C*LC and L = D*LD, where C G W^""^, D G K'^^^''. 

Proof. Since is an interior point, Po ^ and Im(Pi) C Im(Po) for i = 1, . . . , n (See 
[2, Proof of Proposition 2.1].). We have Pq = V*DV, where D is diagonal and V or- 
thogonal. Further on V*LV = i|im(Po) ©O^.^j, d = dim(Im(Po)). Hence L|i,„(pj,) — 
J*(y*LV)J with J* := [/j O'*^^^'*-'^^)] G R^'"'^. Defining Pa Po|im(Po) ^ 0, gives 
Po = B*B, where Po,P G M.'^'^^ and B is invertible. So (p-i)*L|i„,(p,)P"i = 
{B-')*B*BB-' + j:^{B-'rP,U^P^,-^B-'x, - / + E.(S"')*^.|im(p„)Pf 'a;^ L. 
L is in R^""^ and Dl{1) = Di{l). With C* := {B-^)*J*V* G M^xd^dxd^dxd ^ 
R^""^ and D* := (y-^)* JP* G M''^'']R'^^%'^~^'^^ = R'^''^, the lemma is proved. □ 

Theorem 3. Suppose Li G §]R'^^''[x] is a diagonal linear polynomial and the 
set Z3ij(l) is hounded. Then for every linear polynomial L2 G SIR^^^[a;] with 



L2\dl-^{i) ^ 0; there are matrix polynomials Aj 
satisfying 



[x] and Bk G 



pdxe I 



k 



BlLiBk- 



Proof. If Z^Li (1) = 0, then by Theorem[l] —1 is in the convex cone generated by the 
diagonal entries of Li. Now L2 can be expressed in the desired form, since the qua- 
dratic module generated by —1 in any ring with involution consists of all symmetric 
elements by the identity 4a = (a + 1)^ — (a — 1)^. If Dl^ (1) = {a}, then we can use 
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Proposition [T] and we are done. If dim Dli{1) = n, then by Lemma 1 WLOG Li 
is monic and Theorem [2] is used. Otherwise we have 1 < dimI?Li(l) =: k < n — I. 
Since D^^ (1) is convex, it hes in some afBne subspace of dimension k. With transla- 
tion WLOG € Dlj(1) and hence the affine subspace is actuahy a vector subspace 
of dimension k. Let B = {e'j^, ej, . . . , e'j,} be the basis of this subspace. B can 
be completed to the basis of R", i.e., B' = {e'^, . . . , e'^,, ej.^;^, . . . , e^}. Standard 
basis {ei, 62, . . . , e„} of M" can be uniquelly expressed by B' and vice versa, i.e.. 



En (i) I 

.7 = 1 '^j ^3 



and e- = YTj=i l^^j^^J^ ^'^^ unique a"-' ,j3 



ducing new unknows as x'^ = X)J=i l^j^j gives also x. 
expressed Xi-s into Li{xi, . . . ,x„), we get Li{x'i, . . . ,x'^). The map $ 



G M. Therefore intro- 



Putting 



. , ^"^j^ /3j"''aj), is bijective and 



defined by $ : (ai,...,a„) ^ (Z^iPfa 
Li{{ai, . . . ,a„)) = Li($(ai, . . . ,a„)). Hence = Di^{l). So and 

Z)£^(l) are in bijective correspondence. Similarly for 15^2(1) and Z?£^(l). There- 
fore Dli{1) C £>l2(1) ^Li(l) - ^La^-*^)- Froni the construction of basis B' , 
x' e £'£^(1) is of the form (x'^, . . . , a;J., 0, . . . , 0). 
Let us write 




Pi X, 




P'ix'i = Li^i{x[, . 



We notice, that Li is still diagonal. For 
L = Li,i(x;,...,x'fe)® [ ''^+1 



©•■•© 


" x' 










~x' 



(which is obviously diagonal), Di_^{l) = Di{l), and with the use of Theorem [T] 
on diagonal entries of Li and each diagonal entry of L, L can be expressed as 
J2j ^j^j + Sfe Bk^iBk- Hence it suffices to prove the statement of the theorem 
for the pair L, L2. 

Analogously as for Li we write L2 as L2 = (^Rq + P-'t^'ij + Yl7=k+i ^'i^'i = 

L2,i{x'i, . . . ,x'j,) + L2a{x'k+i, ■ ■ ■ ,x'n). We have L2,i|d^^ ^(1) >- 0. Since there exists 
an interior point in Z)^^ Lemma [T] allows us to regard Li i as monic. Finally 
Theorem[5]is used for the pair Li i,L2.i- 

It remains to express L2.2(x'f.^i, . . . = R'k+ix'^,+i + 

cording to Proposition[2l Rix'j^ can be expressed with 

L. Since i was arbitrary, we are done. 
To conclude, we got the expression 

L2K, ...,x'J^J2 ^J^i + E BlL,ix[, x'jBk 









-x' 



R'rix'^ with L. Ac- 
Hence also with 



Using a;,- = X]?=i P^j^^^j^ finally get 



L2{xi 



E 



A*Aj + BlLi{xi, . . . , Xn)Bk 



□ 
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4. Strict positivity of L2\dl{i) 

The next thing to be studied is the necessity of positive definiteness in Theorem 
[5J i.e., whether semidefiniteness suffices. We separately study the one- variable case 
from the general diagonal case. 

4.1. One- variable case. 

Theorem 4. Suppose Li{x) = Pq + Pix G §R'^^''[a;] is a linear polynomial and 
the set _Dij(f) has an interior point. Then for every linear polynomial L2{x) — 
Rq+Rix G SM'^^'^[x] with L2\dj^-^{i) ci 0, there are matrix polynomials Aj G M'^^'^ [x] 
and Bk G R''^^ [x] satisfying 

3 k 

Proof. Case 1: If 0^-^(1) is bounded and there exists an interior point in 1)^^(1), 
then according to Lemma[l] we can assume Li and L2 are monic, i.e., Pq = I, Rq = 
I. Since we have just one variable, we may interpret both Li and L2 as NC 
polynomials, or precisely linear pencils. We will first show, that 
-Pli(I) Q £'l2(1) DL^ C ■■ Let us take X G D^^, X G SM"^™, which means 
Li {X) = / (g) + Pi (g) X ^ 0. Or equivallently (g) / + X (g Pi ^ 0. We have to 
show that X G D^^. Since X is symmetric, it can be real ortogonally diagonalized, 
i.e., UXU'^ — D, where U is an orthogonal matrix of size m. After multiplying with 
invertible matrix C/g)/ we get (C/g)/)(/„g)/-|-Xg)Pi)(C/g)/)^ = Im®I+D®Pi >z 0. 
Hence X G -Dli ^ -D G • Now g) / + _D (g Pi ^ is a block-diagonal matrix 
with the blocks of the form I+d^Pi. It follows Ira®I+D®Pi ^ <^4> I+diPi >z for 
i = 1,2, ... ,m di £ -Dli (1) for i = 1,2, ... ,m. But according to the assumption 
di G Dl^W for alH = 1,2, ... ,m and hence X G Dl^- 

To be able to use LP-satz (i.e., Theorem [5] with L2\dl (1) ^ ^ replaced by 
L2\dl-^ ^ and Aj,Bk constant); cf. [2, Corollary 3.7]; for the pair Li,L2, -Dli 
must be bounded. But by [2, Proposition 2.4] this is equivalent to 13^^(1) being 
bounded. So by LP-satz there exist Bk, such that L2 = J2k B^LiBk. 
Case 2: If ZJ^i (1) is unbounded, then it is an interval of the form [a, 00), (—00, a], 
(—00, 00), a G M. With translation we may assume a = 0. 

First we study the case -Dl^(I) = [0, cxd). Since G Dli{1), we have Pq ^ 0. We 
can also show, that Pi >^ 0. To explain: u*Li{x)u = u*Pqu -\- u*{Pix)u = u*Pou + 
xu*Piu. In the case that u*Piu ^ 0, we have x\u*Piu\ > |u*Pou| for x great 
enough. Therefore, if there exists u, such that u*Piu < 0, then lim x ^ Dl-^{1). 

Contradiction. 

Since Pq and Pi are positive semidefinite, we can use Newcomb's theorem [5, 
Theorem 20.2.2] (It is actually made for complex matrices but with a slight mod- 
ification of the proof it holds for real as well.) to simultaneously diagonalize them 
with invertible S, i.e., S*PoS, S*PiS are both diagonal. So WLOG Li is diagonal. 
Analogously for L2. Now we just use Theorem 1 on diagonal entries of Li and each 
diagonal entry of L2 and we are done. 

In the case ^^^(l) = (— oo,0], we have again Pq ^ 0. As above we show Pi ^ 0. 
Since Pq, Pi are semidefinite, Newcomb's theorem [5, Theorem 20.2.2] can be used 
and we proceed as above. 
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In the case _Dlj(1) = (—00,00), we have Pq h and it is easy to show, that 
Pi = 0. Therefore Li{x) = Pq and analogously L2{x) — Rq, where Bq >z 0. Hence 
L2{x) = C*C. □ 

The following example shows, that 1)^^(1) must have an interior point in Theo- 
rem U) 

Example 2. For the non-monic, non-diagonal polynomial Li{x) ~ [ ]. q ], we 
have that the set Dl-^{1) = {0}. Therefore, Dl-^{1) is non-empty and bounded. It 
also holds that the polynomial L2{x) — x is non-negative on Di^-^{1), but there do 
not exist matrix polynomials Aj e M.[x], G R^^"'"[a;], such that L2 = A*Aj + 
Yjk BlLiBk. 

Proof. Since det(_Li) = —x^, Dl^{1) = {0}. It is obvious, that 2721^^,^(1) > 0. 
The proof will be by contradiction. Let us say there exist Aj, Bk, such that 
J2j A*jAj + Y,k BlLiBk = X. Let Bk be of the form [bf\b^^^Y', where bf^ e R[x] 
and Aj e Ili[a:;]. Comparing the expression J^j AjAj + J^k B'^LiBk with x: 



j k 



The coefficient at 1 on LHS equals X]j ^^,0 + Sfe(^i*^o)^' where Aj,o denotes the 

(k) (k) 

free monomial in Aj and bl q the free monomial in . Since on RHS it is 0, 

Aj,Q — b[''l = for all j, k. But then the coefficient at x on LHS is 0, while on RHS 
1. Contradiction. □ 

4.2. General diagonal case. The aim of this subsection is to prove that for very 
special diagonal Li, we can replace the condition L2\dl^(i) >- in Theorem [2] with 
the weaker condition ^2 1/5^^(1) ^ 0. More precisely, we have the following theorem. 

Theorem 5. Suppose the polynomial Li = Po +X]r=i Pi^i ^ SM''^'*[a;] is diagonal 
and the set 1)^^(1) is an n-simplex. If the polynomial L2 = Ra + Y^^=i ^i^i ^ 
§R^^^[a;], satistfies the condition i2|_Di^(i) b 0, then there exist matrix polynomials 
Aj G R^^^,Bk G R''-^^, such that the following is true: 



3 

L2=Y,A*A,+Y,BlLiBk 

3 k 

In the proof, we will need the following well known result. 

Theorem 6 (1, Theorem 4). Let A be commutative C* -algebra and t : A ^ S a 
positive linear function, where S is a vector subspace of the algebra of all bounded 
operators on some Hilbert space. Then r is completely positive. 

Proof of Theorem\B[ An n-simplex in M" is an intersection of n -I- 1 halfspaces. 
Therefore, it can be defined as Dl{1) of L — 0"^^^ (^a^^ -I- '^j'^^i) = Po + 

ELi^«^' ^ §M("+i)^("+i)[a;], for appropriate af G M. By Theorem 1, L = 
A^LiAfc. Hence, it suffices to prove the statement for the pair L,L2. 
There is an interior point v — {vi,V2, ■ . ■ ,v„) G M" in the n-simplex 
With substitutions the interior point v of Dij(l) becomes the in- 

terior point for L = Pq + J27=iBi^i- ^ ^1^° diagonal, Pq — Lii'^) and 
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Di{l) = Dl{1) - V. The same is true for L2 and Z2. Since Dl{1) C <^ 
Dl{1) - V C Dl,{1) - V ^ Diil) C Di^{l), we have i2bi(i) >= 0. Since Pq ^ 

Li{v), Po is invertible and we have (^(Pq)"') Zi(Po)"^ ==((^'o)"^) Pa{Par^ + 

E^ {(Poy^y iP^S:^){Po)-^ = I + E^ =: Hi), whcrc D^{1) = D^{1). Since 
is an interior point in D^_^{1), by Lemma [1] there exists monic L2, such that 
Di^{l) = Di^{l) and L2 = C*L2C\ L2 = D*L2D. Therefore, it sufHces to prove 
the statement for Li,L2. 

Now we define vector spaces Si := Lin |/, Pi, . . . , P„| and ^2 := Lin |/, Pi, 

...,P„|. Since Dj^{l) is bounded, |/, Pi, . . . , P„| is hneary independent by [2, 

Proposition 2.6] in DR("+^' = {diagonal matrices of size n + 1}. Hence, it is also 
its basis, which implies 5*1 is algebra. Therefore, t : Si S2, where I 1-^ I and 
Pi I— ?> Ri, is a well-defined unital linear map. By [2, Theorem 3.5], it is also positive. 

Since all matrices in |/, Pi, . . . , P„| are diagonal, ^i is commutative algebra. Let 
Si be complex linear span of |/, Pi, ■ • • , Similarly for 5*2, S*^. Now we extend 

T to : S'l ^ §2, where T^{I/Pi) = I/Ri. Since positive elements from Si are 
in 5*1 and T^'j^j^ = r, is positive. Taking Si as A and Sf as S in Theorem 
ini t"' is in fact completely positive. Also, t^I^^ t is completely positive. By 
[2, Theorem 3.5] Dl, C Dl^. By LP-satz [2, Corollary 3.7] for the pair Pi,L2, 
there exist Vj G R'^''"' and fJ. e N, such that L2 = ^^=1 ^/-^i^j- The theorem is 
proved. □ 



Remark. The matrix polynomials Aj, Bk in Theorem [5] are constant while in The- 
orem [2] they are matrix polynomials. 

Remark. If we replace the expression the set Dl^{1) is an n-simplex in Theo- 
rem [5] with the expression the set Lin {Pq, Pi, . . . , P„} is an algebra and has an 
interior point, the theorem still holds. Indeed, since the set Lin{Po,Pi, . . . ,Pn} 
is algebra, it is isomorphic to a subalgebra A in K"+-'^. Further on, we may as- 
sume that this subalgebra separates the n + 1 components of the vector (i.e., given 
1 < i < j < n + 1, we find an element of this algebra, that has distinct i-th and 
j-th component). This assumption does not harm, since otherwise there is redun- 
dancy in Li with respect to the polyhedron defined by Li. By identifying M"+^ 
with the set of all continuous functions from the set X — {1,2, ... ,n} to K, i.e., 
with C{X), and using Stone Weierstrass theorem, we conclude that the algebra 
A is dense in M""*"^. But because of finite dimensionality, it is then equal to the 
full algebra M""*"^. Therefore, the set {Po,Pi, . . . ,P„} is lineary independent and 
since the polyhedron defined by its elements has an interior point, it is an n-simplex. 
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